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We study the total cross section and angular distribution in Rayleigh scattering by hydrogen 
atom in the ground state, within the framework of Dirac relativistic equation and second-order 
perturbation theory. The relativistic states used for the calculations are obtained by making use 
of the finite basis set method and expressed in terms of B-splines and B-polynomials. We pay 
particular attention to the effects that arise from higher (non-dipole) terms in the expansion of the 
electron-photon interaction. It is shown that the angular distribution of scattered photons, while it 
is symmetric with respect to the scattering angle 9=90° within the electric dipole approximation, 
becomes asymmetric when higher multipoles are taken into account. The analytical expression of 
the angular distribution is parametrized in terms of Legendre polynomials. Detailed calculations 
are performed for photons in the energy range 0.5 to 10 keV. When possible, results are compared 
with previous calculations. 
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I. INTRODUCTION 

Rayleigh scattering (also called coherent scattering) de- 
notes the elastic scattering of photons by bound electrons 
and is usually mentioned in order to explain the blue sky 
and the red sunset [l|, [2[ . 

Rayleigh scattering, apart from being a subject of the- 
oretical study by itself, has also interesting applications 
to different fields. Recently, with the availability of X- 
Ray polarization detectors as well as synchrotron and 
FEL sources, new experimental information is emerging 
in the X-ray regime, and, consequently, the demand for 
accurate theoretical predictions is arising [j-Q. 

The total cross section for Rayleigh scattering has been 
widely investigated during the past decades, both within 
relativistic and non-relativistic frameworks [61-flOj. For 
instance, in 2007 Nganso and Njock have presented an- 
alytical results for Raman and Rayleigh scattering in 
hydrogenlike ions based on fully relativistic wavefunc- 
tions and the Sturmian expansions of the Dirac-Coulomb 
Green function [ll| . More recently, Costescu et al. have 
analyzed Rayleigh scattering amplitudes in ions and neu- 
tral atoms by using the independent particle model [12] ■ 

In contrast to the total cross section, less attention has 
been paid to the angular distribution of the scattered 
photons (i.e. the angle-differential cross section) as well 
as their polarization properties. To the best of our knowl- 
edge, a fully relativistic treatment of the angular distri- 
bution in Rayleigh scattering by hydrogen atom based on 
Dirac equation has not been performed yet |13l - [l7| . Hy- 
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drogen atom, although is not, perhaps, the best choice 
for experimental investigation, is the most abundant el- 
ement in the Universe and therefore is of great interest 
in astrophysics. In the extended atmosphere around a 
giant star, for example, near-ultraviolet photons can be 
significantly scattered by atomic hydrogen [18 |. 

In the present work, we study the total cross section 
and the angular distribution of the scattered photons in 
Rayleigh scattering by atomic hydrogen in the ground 
state. The calculations are carried out within a rela- 
tivistic framework through the use of finite basis sets 
for the Dirac equation constructed from B-splines and 
B-polynomials. The angular distribution is furthermore 
parametrized in terms of Legendre polynomials and the 
resulting coefficients are plotted against the photon en- 
ergy. The photon energy range we investigate is 0.5 to 
10 keV. By comparing with previous works, good agree- 
ment is found in analyzing the total cross section in the 
whole energy range. We show that higher multipoles in 
the expansion of the electron-photon interaction opera- 
tor play an important role both in the total cross section 
and, especially, in the angular distribution of the scat- 
tered photons. While, for the total cross section, the non- 
relativistic dipole approximation is adequate for energies 
below ~ 3 keV, in the angular distribution non-dipole 
effects become important already for photon energies > 
500 eV. 

This article is structured as follows. In section [ill we 
present the geometry we consider for the scattering pro- 
cess and the notation used. In section Hill the essential 
theoretical background needed for the calculations is pre- 
sented together with a detailed explanation of the theo- 
retical quantities that we intend to investigate, namely 
the total cross section and the angular distribution. In 



section IIVI wc explain the numerical method we use in 
order to obtain the total cross section and the angular 
distribution. In section |Vl we present our results, by 
showing the impact that high multipoles in the expan- 
sion of the electron-photon interaction have. Finally, a 
summary is given in section IVll together with a few per- 
spectives for further theoretical studies. 



II. ATOMIC SYSTEM AND GEOMETRY 

To explore the cross section and the angular distri- 
bution in Raylcigh scattering, let us first introduce the 
atomic system and the geometry under which the distri- 
bution of scattered photons is considered. 
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FIG. 1: (Color online) The adopted geometry for the scat- 
tering process is displayed. The scattering angle 6 uniquely 
defines the direction of the outgoing photon in the scattering 
[xz) plane. The hydrogen atom is placed at the origin of the 
coordinate axes. 

Our system consists of a hydrogen atom in the ground 
state which is irradiated by light, as shown in Fig. [T] In- 
cident (scattered) light has energy -E-y^a) — ^1(2) j prop- 
agation vector ki(2) and polarization vector £1(2), where 
h is the Planck constant and lo is the angular frequency 
of light. Here and in the following, A denotes the unit 
vector A/|A|, for any vector A. 

We adopt the quantization (z) axis along the direction of 
the incident light (ki). As we will see in Sec. IIIIBl such 
a choice of quantization axis simplifies the multipole ex- 
pansion of the electron-photon interaction operator. The 
scattered light propagates along the direction k2 at angle 
9 with respect to the z axis. The scattering plane {xz) is 
defined by the incoming and outgoing photon directions. 

Now that we have explained the adopted geometry and 
notation, we are ready to present the theory of Rayleigh 
scattering. 



tering both in non-relativistic and relativistic regimes 
j6l-[l]|. In those studies, the analysis of the scattering 
cross section is usually traced back to the second-order 
bound-bound transition amplitude. The procedures to 
derive such an amplitude are explained in standard text- 
books of Quantum Mechanics [l^, [201 , and are based on 
second-order time-dependent perturbation theory. With 
some straightforward algebraic manipulation, the transi- 
tion amplitude for Rayleigh (and Raman) scattering can 
be cast in the form 
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where Wj^i = {Ei, — Ei)/h is the transition frequency be- 
tween states {v) and \i), oc is the vector of Dirac ma- 
trices and E^u\ is the energy of the intermediate (ini- 
tial) atomic bound state. Here, the transition opera- 
tor £1.26*''^'^'" describes the relativistic electron-photon 
interaction. As indicated in Eq. ([T|), the summation 
over the intermediate states runs over the complete one- 
particle spectrum |i/), including a summation over the 
discrete part of the spectrum as well as the integration 
over the positive and negative energy continua. The ini- 
tial state \i) and final state |/) of atomic hydrogen have 
well-defined angular momentum j, angular momentum 
projection rrij and parity (— 1) , where I is the orbital 
angular momentum of the larger component of the Dirac 
spinor. In the following, we will rewrite them respectively 
as I Pi , ji , TTij. ) and | /3y , jj , ruj^ ) , where /3 is a collective la- 
bel used to denote all the additional quantum numbers 
needed to specify the atomic states but for j and m.j . In 
particular, /3 refers to n (principal quantum number) and 
I. 



Absorption First 

S / 




Emission First 
ki.ei k2,e2 




\1\ \i\ 




FIG. 2: (Color online) Feynman diagrams which correspond 
to the first and second term of the transition amplitude ([1} , 
respectively. 



III. THEORY 

A. Evaluation of the transition amplitude 

There has been a lot of theoretical studies on the cal- 
culation of the transition amplitude in Rayleigh scat- 



The transition amplitude in Eq. ([T]) can be interpreted 
in the language of Feynman diagrams from Quantum 
Electrodynamics theory. The whole scattering process 
can be described by the two Feynman diagrams shown 
in Fig. [21 which will be hereinafter called "absorption 
first" and "emission first" , respectively [2l| . The middle 
lines in the diagrams refer to a set of virtual intermediate 



states that represent both electron and positron states. 
Absorption first diagram corresponds to the first term of 
the transition amplitude ([T]) while emission first diagram 
corresponds to the second term. 

Due to the conservation of energy, the quantities Ej-^ ^ 
and Ef^i are simply related by the equation 
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Since Rayleigh scattering is an elastic process, the initial 
and final states are the same, \i) ~ |/), and, thus, Eq. ([2]) 
simplifies to E^-^ = E^^ = E^. 

For gaining deeper insights in the transition ampli- 
tude ((J), in the next subsection we shall decompose 
the electron-photon interaction operator in terms of its 
spherical tensor components. 



where TJ^(f) is defined as 
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The spin spherical tensor ^^ is defined by 
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B. Multipole decompositions of the photon fields 



Since we want to study the angular properties of the 
scattered photons, we need to rewrite the vector plane 
wave ee' "" in terms of elements with well defined angu- 
lar momentum properties. This can be done by using 
the multipole decomposition of the vector plane wave in 
terms of spherical tensors [22|. Such a decomposition 
reads 
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xDi,^iipk,9k,0), (3) 



As seen from Eq. ([3]), the angular dependence of the 
photon emission results from the Wigner (rotation) ma- 
trices. The Wigner matrices transform the multipole 
fields, with original quantization axis along the photon 
propagation direction, into the fields with quantization 
axis along the z||ki direction. Such a choice of quanti- 
zation axis allows us to describe the second photon di- 
rection by means of a single polar angle 6. On account 
of this, we will consider Z?,f„^^(k2 ^ z) = dj^j^-^ (6) and 
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amplitude ([T]). 



where 



aL/(fc'i") = 



AS(fc,r) ifp = 0. 






(4) 



(fc,r) ifp= 1 , 



with A being the photon hclicity. We have defined 
[Li,i2,...,i„] = (2Li + l)(2L2 + l)...(2L„ + l) for the 
sake of simplicity. Each term a'^^j(fc,r) has angular mo- 
mentum L, angular momentum projection M and parity 
(—1)^+^+^. The standard notation A^^Y" is used to de- 
note the electric (e) and magnetic (to,) multipole fields. 
Each of these multipoles can be expressed in terms of the 
spherical Bessel functions jhikr) and the vector spherical 
harmonics T j^ (f ) of rank L as 

Ai™J(fc,r) = j-,.(fcr)Tii(f) , 
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In the following subsection, by substituting the multi- 
pole decomposition ([3]) into the transition amplitude ([Ij , 
we shall show how the total cross section and the angu- 
lar distribution can be obtained in terms of the transition 
amplitude. 



C. Total cross section and angular distribution 



The general form of the relativistic transition ampli- 
tude for Rayleigh scattering is described by Eq. ([IJ . Such 
an amplitude contains the whole information on the prop- 
erties of the scattered radiation. For instance, the angu- 
lar distribution function can be written in terms of the 
squared transition amplitude as (SI units) 
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where a (a = -— ~ ) is the electromagnetic couphng constant and c is the speed of hght in vacuum. Since 
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we do not investigate polarizations in the present article, in Eq. ([8]) we have summed over the final and averaged over 
the initial atom and photon polarizations. By combining Eqs. ([T|), ([3|) and ([H]) and by employing the Wigner-Eckart 
theorem [23|, the transition amplitude can be written as 



Li,L2 PlP2 

Ah 

X E(-^)'-'"'^ (2^-^ + 1)1/2 (e>(l,2)^>(l,2) + e>(2,l)5>(2, 1) 
The reduced (second-order) matrix clement is given by 
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and S-'"{2, 1) is obtained from ((TU)) by i) interchanging the label 1 with 2 and ii) replacing the positive sign in the 
denominator with a negative sign. This latter replacement is given by the fact that, while the second photon is 
emitted, the first photon is absorbed by the atom. 
Following the notation used in [2J] , we have furthermore defined 

e^"-'(i,2) = V(-i)"^/+'"-(2> + i)i/2f ^f \^ ^^\( ^^ ^^ -^O , (11) 

where Q-'''{2, 1) is obtained from Eq. ([TT]) by replacing Li o L2 and Ai -H' M2- 

The transition amplitude ([1]) formally includes the infinite summation over all the multipoles combinations 
PiLiP2L2, where p = 1 (0) refers to E (M). Using standard notation, we may write 

_A/f-77(i ^ /) ~ £:i^i + ElMl + MlMl + MlEl + E2EI + ... . (12) 

However, according to parity and angular momentum selection rules, the number of allowed (i.e. non-zero) terms in 
the summation is restricted to certain combinations of the indices Li, L2, pi and p2. Since the atomic (initial and 
final) state that we consider for the calculations is lsi/2 (the ground state), the allowed multipolc terms are the ones 
with equal parity and \Li — L2I < 1- More explicitly, the multipoles that must be considered are ElEl, A/lAfl, 
£;iM2, M2E1, ... . 
Finally, integrating (|S]) over the second photon directions, we get the total cross section as (SI units) 
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Comparing Eqs. (J13p and (|S]), we note that the sum- Af2), while is inside the modulus squared in the angu- 
mation over the multipoles of the outgoing photon (L2 , 



lar distribution, comes outside the modulus squared in 
the total cross section. By using the notation intro- 
duced in Eq. ([T2|) . we may better underline the con- 
sequences of this fact: Interference terms of the type 
(^1^1)(M1M1)* or {MIM1){E1E1)* are forbidden in 
the calculation of the total cross section, while they 
are allowed in the calculation of the angular distribu- 
tion. Thus, higher (non-dipole) multipoles will start giv- 
ing non-zero contributions in the angular distribution at 
lower photon energies than in the total cross section. 
This effect will be analyzed in detail in Sec. |Vl where 
expressions ^ and P^ will be used to explore the total 
cross section as well as the angular distribution in the 
Rayleigh scattering by hydrogen atom. However, prior 
to doing that, we shall explain how the reduced matrix 
elements ([TU]) are calculated in the present work. 



is ribp - 1). 

The parameters used in both basis sets were optimized 
in order to obtain stability and agreement of six dig- 
its between the results of both basis sets. The opti- 
mal parameters are; i?bs = 60 a.u., nbs = 60, k = 9, 
Rhp = 50 a.u. and Tibp = 40. Such set of parameters was 
already obtained for the case of two photon emission in 
Rcfs. [37l.l38|. 

The finite basis set method has been widely used dur- 
ing the past years to explore the two-photon decay of 
the metastablc 2si/2 state in heavy hydrogcnlike ions 
[21, [3^. In contrast to the two photon decay process, 
to our knowledge no one has ever applied this method to 
calculate the transition rate for Rayleigh scattering. 



IV. COMPUTATION 

During the last decade, various methods have been in- 
vestigated for calculating the reduced second order am- 
plitude ffD|) as well as the transition amplitude ([T]) [lll - 
Il3l . I25I - I271 I . In practice, of course, the summation over 
the complete spectrum contained in ([TU]) is difficult to 
be performed explicitly. Several approaches and approx- 
imation techniques have been proposed to perform such a 
summation. Among them, the Coulomb-Green function 
approach has been widely used for investigating both de- 
cay of and scattering by atoms and ions [2^ (23] . 

An alternative approach is the finite basis set method 
[30l - [35| . The finite basis set method is based on the sup- 
position that the ion (or atom) is enclosed in a finite 
cavity with a radius R large enough to get a good ap- 
proximation for the wave functions. Such a restriction 
leads to a "discretized" continuum part of the atomic 
or ionic spectrum, and hence to a representation of the 
Dirac wavefunctions in terms of pseudo basis set func- 
tions. This basis set forms a complete set of orthonormal 
functions 
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In the present work, we calculate the transition ampli- 
tude (|10p by using B-splines and B-polynomials as finite 
basis sets. The B-splines are one of the most commonly 
used family of piecewise polynomials, since they are well 
adapted to numerical tasks 13 ill . The B-polynomials, or 
the Bernestein polynomials |32| , are a good alternative to 
the B-splines since they allow for analytical finite basis- 
set calculations. These are polynomial functions of nth 
degree that are used to obtain the solution of some linear 
and nonlinear differential equations [32|. The details of 
these basis sets, as well as a comparison between them, 
can be found in Ref. [331 ■ Thus, we restrict ourselves to 
describe the characteristic parameters used in this work. 
The parameters of the B-splines basis set are the radius 
of the cavity (i?bs), the number of B-splines (nbs) and 
their degree (k). As for the B-polynomials, the param- 
eters are the radius of the cavity (i?bp) and the number 
of B-polynomials (ribp) (the degree of the B-polynomials 



V. RESULTS AND DISCUSSION 

With the formalism developed above, we can now ad- 
equately analyze the total cross section and the angu- 
lar distribution in Rayleigh scattering by hydrogen atom. 
The initial and final state considered for the calculations 

is lSi/2. 

In Fig. [31 we plot the total cross section in the photon 
energy range 0.5 to 10 keV as given by Eq. P^ . Calcu- 
lations obtained within the electric dipole approximation 
(ElEl) and by including all multipoles are separately 
displayed. Furthermore, results from Veigele's work |39j 
have been interpolated and are also shown. Our cross sec- 



Veigele 

all multipoles 

E1E1 




E (keV) 



FIG. 3: (Color online) Total cross section in Rayleigh scatter- 
ing by atomic hydrogen, as a function of the photon energy. 
Calculations obtained within the electric dipole approxima- 
tion (red-dashed line) are compared with those including all 
of the allowed multipoles (black-solid line) and with those 
from Veigele [33] (blue-dot-dashed line). 
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FIG. 4: (Color online) Angular distribution in Rayleigh scattering by atomic hydrogen in the ground state as a function of the 
scattering angle 0. Results are calculated within the exact relativistic theory (solid-black line) and within the electric dipole 
approximation (red-dashed line), for three selected photon energies. 



tion practically coincides with Veigele's, indicating that 
there is good agreement between the two calculations. 
Moreover, it can be seen that for high photon energies 
the electric dipole approximation deviates from the re- 
sults of Vcigele and of us obtained with the account of 
all multipoles. For photon energy w 3 kcV, the electric 
dipole approximation underestimates the cross section by 
about 3%, while, for photon energy > 6 kcV, it is lower 
by a factor of two or more. Indeed, this result is not un- 
expected: in low photon energy range, the first (dipole) 
term ElEl dominates the transition amplitude, while, 
for higher photon energy (such that fcr > 1, where r is 
the atomic radius), higher multipoles play an important 
role, as it is evident from Eqs. ([3]) and ([5]). For a de- 
tailed comparison, our results for the total cross section 
are presented in Table |T] together with Veigele's [33 and 
NIST's (National Institute of Standards and Technology) 
[401, for selected photon energy values. 

In addition to the total cross section, we analyze the 
angular distribution of the scattered photons. In Fig. SI 
wc display the angular distribution as obtained from 
Eq. ([8]), for the three photon energy values E^= 500 cV, 2 
kcV and 5 kcV. Calculations have been performed within 



TABLE I: Different calculations for the total cross section for 
Rayleigh scattering by hydrogen in the ground state, for se- 
lected values of the photon energy. The results are compared 
with Veigele ^ and NIST [13] . 



E., (eV) 


Veigele 


NIST 


This calculation 


500 


0.642 


n.a. 


0.6465 


1000 


0.579 


0.5805 


0.5812 


2000 


0.414 


0.4141 


0.4142 


3000 


0.277 


0.2760 


0.2764 


5000 


0.135 


0.1341 


0.1341 


8000 


0.0618 


0.0612 


0.0609 


10000 


0.0416 


0.04121 


0.04030 



the electric dipole approximation and by taking into ac- 
count all multipoles that have a contribution greater than 
1 percent. As seen from the figure, the angular distribu- 
tion for low photon energy values is well described by 
the clcctric-dipole approximation and by the expression 
~ 1 -f- cos^ [a| , as expected from the non-rclativistic 
theory. For photon energy values w 2 kcV, forward scat- 
tering becomes dominant, compared to backward scatter- 
ing, and leads to an asymmetric angular distribution. For 
even higher photon energy values (> 5 keV), the back- 
ward scattering is strongly suppressed and « 90% of scat- 
tered photons are found within the interval < 6 < 60°. 
This effect is known from previous calculations and ex- 
periments on other ions and atoms. See for instance 
Rcf. 0- 

By comparing Figs. [3] and 21 we notice that the an- 
gular distribution is much more sensitive to non-dipole 
effects than the total cross section: While, for the total 
cross section, non-dipole effects are suppressed below ~ 
3 keV, in the angular distribution they become already 
important for photon energies > 500 eV. As discussed 
previously (see Eq. ([T?)) and the paragraph afterwards), 
this is a direct consequence of the different structure 
of the amplitudes M and M displayed in Eqs. © and 
([T4)) . respectively. In particular, following the notation 
used in Eq. (J12p . the main contribution for the asym- 
metric shape of the angular distribution comes from the 
(£:i^l)(A/lA/l)* and (A/lil/l)(£'l£'l)* terms, which 
are forbidden in the calculation of the total cross section. 

Figure m displays the angular distribution for just three 
specific energy values. To present the angular distribu- 
tion in the whole energy range, we make an expansion of 
dc/dfi in terms of Legendrc polynomials and wc normal- 
ize it to the total cross section a: 



^ = ^(/3oPo(cos0) +/3iPi(cos0) 



(15) 



where (3i are real numbers called anisotropy coefficients 

The results for the (3 coefficients are displayed in Fig. [SI 
within the considered energy range. The first coefficient 
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least for photon energies far from the ionization thresh- 
old. The angular part of the Dirac wave function remains 
unchanged when the atomic number Z is increased, while 
the radial part gets contracted: The mean radius, r, of 
the electronic ground state gets shrunk and scales like 
r ~ r^/Z , where rp denotes the radius of hydrogen atom 
[42l |. The transition operator in equation ([5]) contains 
both the angular operator f and the radial operator kr. 
Thus, while the action of the former will be the same for 
hydrogenlike ions as for hydrogen atom, the action of the 
latter will scale as kr ~ kr^/Z = {k/Z)ro. Eventually 
therefore, for energies far from the ionization threshold, 
the whole angular distribution will scale in the energy 
domain with a factor of 1/Z. More explicitly, we approx- 
imately obtain the angular distribution of photons which 
scatter by hydrogenlike atoms with atomic number Z by 
taking the angular distribution for hydrogen atom and 
replacing k — > k/Z. 



FIG. 5: (Color online) The anisotropy coefficients /3i of Eq. 
(|15|l are plotted against the photon energy. Results obtained 
within the electric and magnetic dipole approximation (red- 
dashed line) are compared with results obtained by taking 
into account all photon multipoles with Li < 4 and L2 < 4 
(black-solid line), which are dipole, quadrupole, octopole and 
hexadecapole moments of each photon field, both of magnetic 
and electric type. 



(3q is constant for any photon energy and therefore is 
not displayed; /3o ~ 0.0796. The red-dashed curves are 
obtained by including only the electric and magnetic 
dipole contributions (ElEl+MlMl) while the black- 
solid curves are obtained by taking into account the mul- 
tipoles with Li < A and L2 < 4. Those multipoles are 
nothing but dipole, quadrupole, octopole and hexade- 
capole moments of each photon field, both of electric and 
magnetic type. From the figure, one immediately notices 
that, as expected, the lower the photon energy is, the 
fewer coefficients are needed for an adequate description 
of the angular distribution. 

As seen from Fig. [51 for low photon energies we have 
/3o ~ 2/^2 and /3i ss for i 7^ 0, 2. By explicitly writing the 
definitions of Legendre polynomials, one can easily check 
that the well-known shape ^1-1- cos^ 9 for the angu- 
lar distribution is obtained within the (non-relativistic) 
Thomson limit. In contrast, for higher photon energies, 
higher-order Legendre polynomials must be invoked to 
well describe the shape of the angular distribution, as 
the number of multipoles that contribute to it increases. 

Although the calculations have been performed for hy- 
drogen atom, the behavior of the angular distribution 
for hydrogenlike ions can be qualitatively predicted, at 



VI. SUMMARY AND OUTLOOK 



In summary, we studied the total cross section and 
angular distribution in the Rayleigh scattering of pho- 
tons by hydrogen atom, based on Dirac relativistic equa- 
tion and second-order perturbation theory. We decom- 
posed the transition amplitude in terms of spherical ten- 
sors and (radial) reduced amplitudes, where these latter 
include a summation over the whole atomic spectrum. 
This summation was performed by means of finite basis 
set method. The multipole expansion of the photon fields 
allowed us to investigate non-dipole effects both in the to- 
tal cross section and in the angular distribution. As for 
the total cross section, we compared our results with pre- 
vious calculations and good agreement has been found. 
We then studied the angular distribution of photons in 
the energy range 500 eV to 5 keV. We here found that 
non-dipole effects become important at much lower en- 
ergies than for the total cross section. As discussed, this 
has to be attributed to the different structure of the tran- 
sition amplitude for angular distribution and total cross 
section. Finally, in order to give a more practical account 
of the angular distribution function, we expanded it with 
Legendre polynomials and we plotted the resulting co- 
efficients against the photon energy, noticing an evident 
convergence. 

Owing to recent advances in detector technology [43{ . 
photon polarization studies in atomic processes have 
lately become important [44| - |46| . Further studies are in 
progress to investigate the polarization properties of the 
scattered photons. 
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